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Abstract. - The dynamical properties of the finite-size magnetization M in the critical 
region T < Tktb of the planar rotor model on a L x L square lattice are analyzed by means of 
microcanonical simulations . The behavior of the q = structure factor at high frequencies is 
consistent with field-theoretical results, but new additional features occur at lower frequencies. 
The motion of M determines a region of spectral lines and the presence of a central peak, which 
we attribute to phase diffusion. Near Tktb the diffusion constant scales with system size as 
D^L- 1 - 6 ^. 



The XY or "planar rotor" model consists of a set of classical spins S r of unit length confined 
in a plane, whose orientation is specified by the angle r , r being the position vector on a 
square lattice of size N = LxL with periodic boundary conditions. It is known that this model 
does not admit equations of motion and therefore its canonical dynamics is usually simulated 
either by Monte-Carlo methods fl|, || or by Langevin type equations ||, ||. Microcanonical 
approaches consist instead, either in considering a three component spin model || or into 
adding a kinetic energy term The latter method can be also generally applied to other 
physical systems (see e.g. the application to lattice gauge theories 0). All these methods 
should display the same static properties, as it has been verified up to some extent ||§|,p||. 

One of the most striking features of the XY model is the presence of strong finite-size 
effects e.g. the existence of a sizeable magnetization for large samples, despite the fact 
that long-range order cannot occur for the infinite system. For this reason, we are motivated 
to study the dynamical features of a finite lattice, a topic which has been far less investigated. 
In this respect the choice of the dynamics, e.g. canonical vs. microcanonical, is expected to 
play a crucial role. 

In this Letter we consider the following Hamiltonian 

n = Ey+ Et 1 -™ 8 ^'-^)], (i) 
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where p r — 9 r is the angular momentum of the rotor and the sum ranges over the four nearest 



neighbors of site r 1 1 1 . We have set both the inertia of the rotors and the ferromagnetic 
coupling constant to unity so that the only physical control parameter is the energy per spin 
e = 7i/N . Actually, there is a second constant of the motion, the total angular momentum 
P = ^2 r p r , whose choice affects the results in a trivial way. In the numerical simulations we 
set P = to avoid global ballistic rotation. 

A previous study of the static properties of ( |l|) showed that the system undergoes a 



Kosterlitz-Thouless-Berezinskii (KTB) transition |12| at e = cktb ~ 1-0 corresponding to a 
kinetic temperature Tktb ~ 0.89, which is in agreement with the value 0.894(5) obtained in 
the canonical ensemble Q . 

Here, we focus our attention on the statistical behaviour of the istantaneous complex 
magnetization M = M x + iM y , defined as 

M = ±J2S r = ±J2e^ =m^, (2) 

r r 

in the ordered phase e < cktb- In this regime the phase ip is well defined, because m is 
bounded away from zero pUfl . Moreover, (m) is a sizeable quantity and can be regarded as 
a pseudo {i.e. non-extensive) order parameter ((•} denotes the ensemble average). Indeed, it 
vanishes only algebraically with system size, (m) = (2A)~ fcsT / 8Tr in the spin-wave limit . 

The numerical integration of the equations of motion is performed using the fourth-order 
McLahlan-Atela algorithm [Q which is an explicit scheme constructed from a suitable trun- 
cation of the evolution operator that preserves the Hamiltonian structure. One of the major 
merits of symplectic algorithms is that the error on the energy does not increase with the 
length of the run. The chosen time step (0.01-0.05 in our units) ensures that in every simula- 
tion energy fluctuates around the prescribed value with a relative accuracy below 10~ 5 . 

The dynamics is fully characterized by the structure factor 5(q, uj), i.e. the space-time 
Fourier transform of the the spin-spin correlation function (we assume translational invariance) 

S(T,t) = ( e *[«'(*)- fl o(0)]) . ( 3 ) 

Accordingly, the power spectrum of M is proportional to the zero wavenumber (q = 0) 
structure factor 

/>oo 

S(0,u) oc / (M(t)M*(0))e iut dt = (\M(uj)\ 2 ) . (4) 



The dynamical finite-size scaling hypothesis suggests that this quantity should be written 
as (\M(w)\ 2 )/AL z = f(uL), where A is the area below the spectrum ||. According to the 
theory by Nelson and Fisher (NF) [jl5| we expect z = 1 in our case. Our numerical results 
are in good agreement with these predictions, see Fig. [l]a. Moreover, the high-frequency part 
(Lu) S> 1) displays a power-law decay over approximatively one decade. This is a prediction of 
NF theory which, relying on a hydrodynamic description, is expected to hold in this frequency 
range, corresponding to times much shorter than the typical transit time of a wave on the 
lattice. Obviously a cutoff frequency of order unity is imposed by the lattice discreteness. 
Although an accurate estimation of the decay exponent is difficult in such a narrow range, 
the available data are consistent with the predicted law u>~ 3+ri (with rj — ksT/2Tr) at 
least up to the energy where a significant density of vortex pairs appears. This is precisely 
the case at e = 0.992 where we found uj~ 2 - b ^ which is actually significantly different from the 
NF value 3 — ij with rj k, ^(Tktb) = 1/4- Notice that this result is somehow in contrast with 
the theoretical expectations [Q as it corresponds to 77 > 1/4. Similar deviations are indeed 
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observed in the static case and are traced back to the presence of multiplicative logarithmic 
corrections. Unfortunately, even the most accurate numerical result are controversial and 
apparently in contrast with the theory (see |l7j | and references within for a thoughtful account 
on this issue). We thus limit ourselves to observe that, at variance with what is reported in 
the literature we find 77 > 1/4 in the critical region. 

The additional structure present in the spectrum at lower frequencies is beyond the validity 
range of NF theory and is related to the long-time dynamics of M for the finite lattice. Indeed, 
it has been previously observed, in both Monte-Carlo |l3| and microcanonical simulations, 
that M fluctuates in time in the complex plane in a narrow region around the circle of radius 
(to) centered in the origin. A closer inspection to our data, reveals that such a motion has 
two main components. On one hand, m(t) oscillates irregularly around its mean value, with a 
typical frequency fl which is inversely proportional to L. On the other hand, the phase ip(t) 
performs a slower random motion whose displacement rapidly decreases with the energy. For 
instance, already below e ~ 0.3 the phase changes of M are hardly observable on the typical 
size and time scales of our simulations (10 3 — 10 4 time units). 

These features of the motion of M reflect into some distinctive properties of S(0,oj). 
More precisely, the oscillations of m manifest themselves as a very sharp line at 51 (Fig. |l|a), 
whereas the slower phase motion determines a low frequency (to <C f2) component centered 
at uj = (Fig. |l|b). The latter is naturally interpreted as a signature of phase diffusion, 
once we assume -0 to be a Brownian variable with a diffusion time D^ 1 much larger than the 
typical time scale of the motion of m. Accordingly, the spectrum is well fitted by a Lorentzian 
curve C/(D 2 + w 2 ), which allows also to determine D. Furthermore, the relative weight of 
the diffusive component of (|M(w)| 2 ) decreases upon decreasing e, in agreement with the 
qualitative observation that phase motion becomes less and less pronounced. 

To further support this interpretation, we have computed the time evolution of the mean 
squared value a 1 = ((ip — 4>o) 2 ) for several lattice sizes. As shown in Fig. ||a, the data are 
consistent with a diffusive motion, i.e. a 2 — 2Dt, for large enough t. We have also performed 
an independent measure of the diffusion constant based on Green-Kubo formula 



Remarkably, the values computed using the three methods agree extremely well within the 
statistical accuracy (better than 15% in the worst case - see the inset of Fig. ||a). 

In analogy with what happens for the static magnetization we expect the diffusion co- 
efficient to approach zero on increasing the lattice size. Indeed, we found evidence that D 
vanishes algebraically with a nontrivial exponent: a best fit to the data in the inset of Fig. || 
gives D cx L _1 - 6 ( 3 ) in the range 8 < L < 128. Although we did not attempt to perform a 
systematic study, we have observed that a similar scaling holds also at lower energies. For in- 
stance, at e = 0.70 we found a sligthly smaller exponent w 1.5 in the same range of sizes. Thus, 
we cannot rule out the possibility that such exponent depends on temperature in analogy with 
the one of finite-size magnetization. 

The oscillations of m are related to the spin-wave modes with the longest wavelength. 
We thus expect to be proportional to c/L, where c is the renormalized spin- wave velocity 
appearing in the effective dispersion relation 




(5) 



2 a 2 r • 2 1x . 2 Qy 

£ = 4c sin h sm — 

q L 2 2 



(G) 



Neglecting the contribution of vortices, c is only determined by the anharmonic interaction 
among the spin- wave modes, thus leading to the estimate —8c 2 log c = fcgT ||. The data 
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Fig. 1 " Fig. 2 

Fig. 1 - (a) Scaled high-frequency part of the spectrum (\M(lu)\ 2 ) for e = 0.992 and L = 64, 96, 128. 
Each curve is an average of ~ 10 3 runs, each of 5 10 3 time units (sampling time 0.3). The thin dashed 
lines represent the power-laws lo~ 2 due to diffusion and the NF decay ui~ 2 ' 75 . (b) Low-frequency part 
of {\M(uj)\ 2 ) obtained from longer runs (3 10 5 time units) for L = 8, 16, 32 (from top to bottom). The 
spectra are not scaled and are in arbitrary units. Thin dashed lines are Lorentzian fits C/(D 2 +lo 2 ). 



reported in Fig. show that this results accurately describes the dependence of on T at 
least for low enough temperatures. Systematic deviation at larger temperature have of course 
to be expected due to the appearance of vortices. Indeed, it is usually believed that their 
effect is to further diminish the spin-wave velocity, which is in qualitative agreement with our 
data. 
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Fig. 2 - (a) Mean squared displacement a 2 = ((ip — i/jo) 2 } at e = 0.992 for L = 16, 32, 64 (from top 
to bottom). The inset shows the corresponding diffusion coefficient D evaluated from a linear fit 
(circles), from Green-Kubo formula (squares) and from a Lorentzian fit of the low- frequency region 
of (\M((jj)\ 2 ) (see Fig. |I|d). (b) Scaled oscillation frequency LQ. versus the kinetic temperature for 
L — 64, 128 (circles and triangles respectively). The solid line is proportional to c(T) (see text). 
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In order to better illustrate the new effects introduced by the finite lattice, we have evalu- 
ated the correlation (^) in the low-temperature (spin-wave) region to lowest order (Gaussian 
theory). The resulting expression 



is the lattice version of the formula obtained by NF in the continuum |15| . Their results can 
by summarized by saying that for ct 3> |r|, S(r,t) should be independent of r and vanish 
algebraically in time as (ct)^ 11 . We have evaluated numerically formula ([?]) by performing 
the sum over the discrete set of allowed qs for several values of r. As a first check, we have 
found that 5(r, t) is practically independent of r for r <C L. Therefore, we show in Fig. || 
the behavior of S(0,t) for several values of L. We observe a dramatic difference with what 
expected from NF theory: (i) the power law decay is restricted to a small region of short 
times (t <C L) and rather well-defined oscillations appear with typical period of order L at 
larger times; (ii) the correlation function does not vanish asymptotically, as it should be for 
a system with nonzero magnetization. The oscillating behavior of (f^) is clearly due to the 
long-wavelength spin-waves and is therefore connected with the line structure around £1 in 
Fig. ^]a. Hence, this signals again that beyond the limit of validity of NF theory new effects 
due to the finite-size appear. 

To summarize, in this Letter we have studied finite-size effects on the Hamiltonian dynam- 
ics of the magnetization vector M in the 2D XY model. We have investigated the consequences 
of the motion of the modulus and the phase of M on the behavior of the q = structure 
factor. They amount to the appearance of a diffusive central component, which we associate 
to phase motion, and of a peak due to the oscillations of the modulus. The high frequency 
part is basically consistent with Nelson-Fisher theory pj|. Similar properties are observed 
for other types of dynamics for which the shape of 5(q, t) deviates from the one given by 
the Nelson-Fisher theory, exactly for the presence of a central peak. We conjecture that such 
features might be related to our results, and a more systematic comparison would be therefore 
desirable. 

Furthermore, we have found a convincing evidence that the diffusion constant D vanishes 
algebraically with L with an exponent different from z — 1. If this would be confirmed for 
larger lattices, one should conclude that the correlation function does not exhibit a single 
scaling with L z , which in turn would imply a violation of the dynamic scaling hypothesis. 
This should be signalled by a bad superposition of the structure factor below u> ~ D for 
different sizes. A direct test of this effect would however require simulations on much longer 
times. Since scaling violations indeed occur for coarsening phenomena pc| , which are far 
off-equilibrium processes, it would interesting to investigate further this possibility also close 
to equilibrium. 
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Fig. 3 - Finite-size scaling of the correlation function S(0,t), formula (|7j) for fcgT = 0.05 and L — 
128, 256, 512, 1024. The scaling exponent is rj = k B T/2n « 0.0079 . . .. 
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